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ABSTRACT 

We present, for the first time, a statistical study of [N ii] 205 /im line emission for a large sample of 
local luminous infrared galaxies using Herschel Spectral and Photometric Imaging Receiver Fourier 
Transform Spectrometer (SPIRE FTS) data. For our sample of galaxies, we investigate the correlation 
bet'wcen the [N ii] luminosity (Ljnu]) and the total infrared luminosity (ija), as "well as the dependence 
of i[Nii]/-^iR ratio on Ljr, far infrared colors (IRAS /eo/Aoo) a-nd the [O ill] 88 /zm to [N ii] luminosity 
ratio. We find that Lj^u] correlates almost linearly with Lir for non AGN galaxies (all having 
LiR < 10^^ Lq) in our sample, which implies that i[Nii] can serve as a SFR tracer which is particularly 
useful for high redshift galaxies which will be observed with forthcoming submm spectroscopic facilities 
such as the Atacama Large Millimeter/submillimeter Array. Our analysis shows that the deviation 
from the mean L[Nn]-iiR relation correlates with tracers of the ionization parameter, which suggests 
the scatter in this relation is mainly due to the variations in the hardness, and/or ionization parameter, 
of the ambient galactic UV field among the sources in our sample. 
Keywords: galaxies: evolution — galaxies: ISM — galaxies: starburst — infrared: ISM 



1. INTRODUCTION 

The star formation rate (SFR) is one of the fundamen- 
tal physical parameters for characterizing galaxies. Hav- 
ing an effective way to derive SFRs for galaxies spanning 
a large range of look-back times is crucial for our under- 
standing of galaxy evolution leading back to initial condi- 
tions. Currently, SFRs have been inferred using contin- 
uum or spectral line emission from a wide range of wave- 
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lengths (e.g. Calzetti ct al. 2009 and references therein; 
Kcnnicutt & Evans 2012). However, few of these SFR 
indicators can be observed over a large redshift range us- 
ing current facilities. Using the capabilities of the ESA 
Herschel Space Observatory (hereafter Herschel] Pilbratt 
et al. 2010), it is possible to define new SFR indicators in 
the far infrared (FIR) window, which can be followed up 
by the Atacama Large Millimeter/submillimeter Array 
(ALMA) over a large redshift range. 

The [Nil] 205 /xm line, which arises from the ^Pi— >-^Po 
transition (205.197 ^m; hereafter [Nli](l-i'0)) of the 
ground state of the singly ionized nitrogen, is expected 
to be an excellent indicator of SFR. This FIR line emis- 
sion essentially traces all of the warm ionized interstel- 
lar medium (ISM). The ionization potential of nitrogen 
(14.53 eV) is only ^ 0.9 eV larger than that of hydrogen 
(13.6 eV), allowing the [Nil] luminosity to be used as 
an estimate of the flux of ionizing photons (e.g. Bennett 
et al. 1994). The critical density of [Nil] (l-J-O) is only 
44 cm-3 at T = 8000 K (Oberst et al. 2006), and the 
energy level that emits the line corresponds to only 70 
K, so that it can be easily collisionally excited. In addi- 
tion, the [N ii] 205 /xm emission is usually optically thin 
because of its small Einstein coefficient, and suffers much 
less dust extinction than optical and near infrared lines. 
In particular, this line is accessible using ALMA for a 
significant fraction of the 0.6 ^ ^ 16 redshift range 
(bands 3-10j3- For more local galaxies, the [Nil] (1— 5>0) 
line can be observed in the two windows around 200 fim 
(1.25 - 1.4 THz and 1.45 - 1.6 THz; e.g. Paine et al. 
2000; Yang ct al. 2010) from the very driest sites with 
future facilities such as the Sub-mm/THz telescope at 

Redshift intervals corresponding to different ALMA receivers: 
11-6 < ^ < 16.3 (Band 3); 3.0 < z < 10.6 (Bands 4 - 7); 2.0 < z < 
2.7 (Band 8); 1.0 < z < 1.4 (Band 9); 0.6 < 2 < 0.8 (Band 10) 
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the Chinese Antarctic Observatory under development 
(e.g. Cui 2010). 

Early studies have shown that the [Nil] 205 fim line 
is comparatively bright. In most normal star- forming 
galaxies, it is the third to fifth brightest FIR/submm 
line after [Cii] 158 /im, [Oiii] 88 ^m, [Oi] 63 ^m and 
[Nil] 122 fim lines (e.g., Wright et al. 1991; Brauher et 
al. 2008). The luminosity of [Nil] 205 /xm line (L[nii]) 
may be up to lO^"' times the total infrared luminosity 
(Lir). For example, -/j[nii] is about 3 x 10~^Lir and 
5 X IQ-^iiR in the Milky Way (Wright et al. 1991) and 
M82 (Petuchowski et al. 1994), respectively. With such a 
high luminosity, this line has already been detected with 
ALMA for luminous/ultraluminous infrared galaxies at 
high rcdshift (e.g. Nagao et al 2012; z = 4.76). 

The [Nil] 205 fim line is generally inaccessible to 
ground-based facilities, and its rest wavelength is longer 
than the cutoff of the Infrared Space Observatory (ISO) 
Long Wavelength Spectrometer (LWS). Therefore, this 
line was observed in only a handful of extragalactic ob- 
jects using satellite and airborne platforms (e.g. Wright 
et al. 1991; Petuchowski et al. 1994; Lord et al. 1995) 
prior to the advent of Herschel. Here we report our 
first results on the [N ii] 205 /xm line emission for a large 
sample of galaxies, observed with the Fourier-transform 
spectrometer (FTS) of the SPIRE instrument (Griffin et 
al. 2010) on board Herschel. We present the correla- 
tion between L[nii] £^nd Lir, to characterize the [Nil] 
205 /Ltm line as a SFR indicator. We also investigate 
possible dependences of [N ii] 205 fim to IR luminosity 
ratio (L[nii]/-^ir) on Lir, FIR color (/eo/Zioo, where /go 
and /loo are the IRAS 60 /im and 100 /im fluxes, respec- 
tively), and the [Oiii] 88 /im to [Nil] 205 /im emission 
ratio (i [OIII] /-^ [Nil]), to reveal the cause of the scatter in 
this relation. Throughout the paper, we adopt a Hub- 
ble constant of Hq = 70 km s~^ Mpc~^, Hm = 0.3, and 
Qa = 0.7. 



2. SAMPLE, OBSERVATIONS AND DATA REDUCTION 

The sample discussed in this Letter is part of the Her- 
schel open time project Herschel Spectroscopic Survey of 
Warm Molecular Gas in Local Luminous Infrared Galax- 
ies (LIRGs) (PI: N. Lu), which aims primarily at study- 
ing the dense and warm molecular gas properties of 125 
LIRGs (LiR = L(8 - 1000 /im) > IO^Lq, where ijR 
was calculated by using the IRAS four-band fluxes and 
the equation given in Sanders & Mirabel (1996)), which 
comprise a flux limited subset of the Great Observato- 
ries All-Sky LIRGs Survey sample (GOALS; Armus et 
al. 2009). The full program details as well as obser- 
vational data on individual galaxies will be given in a 
future paper upon the expected completion of the pro- 
gram in early 2013. Here we present a subsample of 70 
galaxies, including 63 LIRGs and 7 ULIRGs (iiR > lO^^ 
Lq), all of which are point sources with respect to the 
Herschel SPIRE/FTS beam. The determination of these 
galaxies as relative point sources was based both on their 
observed mid-infrared emission extent (for example, see 
Diaz-Santos et al. 2010, 2011), and on the consistency 
between the observed emission strength seen in the over- 
lapping FTS bands between 316 and 324 /im (beam sizes 
are ~ 20" and ~ 37" for SSW and SLW, respectively, in 
this wavelength range). A companion letter (Lu et al. 
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Figure 1. Typical spectra (corrected to the source rest frame) of 
the [N ii] 205 line for galaxies with a range of /go //lOO colors 
(going redder from top to bottom). In each panel, the blue, red 
and black linos show the observed spectrum, the fitted profile of 
the [Nil] line plus continuum, and the residual, respectively. In 
the left panels, the [Nil] line is fitted using a Sine function with 
a free width parameter. The right panels show the same data but 
the line is fitted using a fixed instrumental Sine function convolved 
with a free width gaussian shape (SCO). It is clear that the SCG 
function can give a bettor fitting result for this line. 

2013) gives a brief program description and some initial 
results on the observed CO line emission. 

The observations were conducted with the SPIRE FTS 
in staring and high spectral resolution mode between 
194-672 /im, yielding a resolution of 0.04 cm~^ (=1.2 
GHz). The data were reduced using a customized ver- 
sion of the standard SPIRE reduction and calibration 
pipeline for point source mode included in HIPE version 
7.3 (Ott 2010). We used an averaged telescope relative 
spectral response function (RSRF) in place of the daily 
RSRF, to perform the flux calibration. During the reduc- 
tion, we did not apply any apodization to the spectra. 
To extract the flnal spectrum, a "background" emission, 
which was obtained by fltting the median spectrum from 
the spatially surrounding detectors, was subtracted from 
the central detector. 

In most eases the [Nil] 205 /im line is the brightest 
line in the SPIRE FTS wavelength range, and almost 
all sources show high signal-to-noise ratio (S/N) in the 
[N ii] line. In Fig. 1 we plot several typical spectra of the 
[Nil] line for galaxies with a range of FIR colors (0.4 < 
/eo/Zioo ^ 1-4), and superimpose a fitted continuum and 
the line profile. From the figure, one can see that the 
fixed Sine function convolved with a free width gaussian 
(SCG) profile gives a better fitting result for the [Nil] 
line since the instrumental resolution of the SPIRE FTS 
in high resolution mode at ^ 210 /im is only about 300 
km and the line width has an infiucncc on the line 
shape. Therefore, we adopted the SCG profiles for the 
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integrated [N ii] line flux measurements except for one 
galaxy (CGCG448-020) where a single Sine profile was 
better due to its narrow line width (< 200 km s~^ for 
low J CO lines; e.g. Baan et al. 2008; Leech et al. 
2010) and relatively low S/N 4) in the line. The 
resulting full width at half maximum of the [N ii] line 
(excluding CGCG448-020) is 140 - 610 km s-\ with a 
median value of ~ 300 km s~^. In most cases, the Icr 
statistical uncertainties in the integrated line flux are less 
than 10%, and the median value is 7%. 

3. RESULTS AND DISCUSSIONS 
3.1. i[Nii]--^iR Correlation 

Our sample galaxies are all (U)LIRGs and 23 out of 
70 show AGN/LINER features (identified using optical, 
mid infrared, or X-ray data). In order to increase the 
size and dynamic range of the sample used to study the 
-^[Niil--^iR relation, we also include 30 unresolved galax- 
ies observed by ISO LWS and compiled by Brauher et 
al. (2008), in our analysis. For these objects, most of 
which are normal star-forming galaxies, the luminosity 
distances were derived with the same method as used 
for our (U)LIRG sample (e.g. Armus et al. 2009), i.e. 
by correcting the heliocentric velocity for the 3-attractor 
flow model of Mould et al. (2000). Their [Nil] 205 //m 
fluxes arc derived from the observed [Nil] 122 /^m fluxes 
given in Brauher et al. (2008), using the theoretical [Nil] 
122 /xm to [Nil] 205 /im emission ratio (hereafter R21) at 
Uf. = 80 cm""^ (the median value for late type galaxies; 
see e.g. Ho et al. 1997) of i?2i=2.6. Here we compute 
the theoretical i?2i by assuming electron impact exci- 
tation for the ground-state levels of N+, and using the 
fitted resuhs (Draine 2011) of Hudson & BcU (2005) with 
T = 8000 K to calculate the collision strengths (more de- 
tails are provided in Zhao et al. 2013, in preparation). 
Regarding the inclusion of Brauher et al. galaxies in our 
study, we note that (1) about 90% of the galaxies in Ho 
et al. (1997) have Ue < 300 cm~'^, yielding a theoretical 
R21 = 0.94-5.3 (rie = 10-300 cm^^); (2) empirically, we 
see an observed i?2i of 1.5— 4.3 for the nine galaxies (eight 
in our sample plus M82) that have both [Nil] 122 /im and 
[Nil] 205 fim data. Therefore, the uncertainty in i[Nii] 
caused by extrapolating [Nil] 205 /im fiuxes from [Nil] 
122 /im emissions is a factor of a few (~ 0.3 dex), and 
is significantly smaller than the two order-of-magnitude 
span of LiR found in our sample, thus yielding useful 
results concerning the Ltnuj-Lir relation. 

In Fig. 2, we plot ivjNn] against Lir for both our 
(U)LIRGs (cirles) and Brauher et al.'s galaxies (dia- 
monds), as well as the starburst galaxy M82 (star- 
symbol; Pctuchowski et al. 1994). The solid symbols 
indicate AGNs/LINERs. Generally, increases with 

LiR, albeit the scatter seems to increase with Lir (we will 
return to this later). In fact, the scatter at the high Lir 
end will be much reduced if we only look at the non- 
active galaxies. A nonweighted least-squares linear fit, 
using a geometrical mean functional relationship (Isobe 
et al. 1990), to the non-active galaxies, gives 

logiiR = (4.51 ± 0.32) + (0.95 ± 0.05) logi[Nii] (1) 

with a scatter of 0.26 dex in Xir. The Spearman rank 
correlation coefficient p of the trend is 0.78 at a > 5a level 



of significance, which indicates a very strong correlation 
between I/[nii] and Lir. Our fit shows that, on average, 
LiR scales with i]^jj] with a between 0.8 and 1.1 at the 
3(7 significance. This nearly linear relation suggests that 
the power source of the [Nil] 205 /im emission tightly 
relates to the star-forming activities. Using the SFR- 
LiR calibration given in Kennicutt & Evans (2012), we 
can also build a direct relation between SFR and i[Nii], 
namely 

logSFR(MQyr-i) = (-5.31±0.32)+(0.95±0.05) logL[Nii] 

(2) 

Comparing to the IR luminosity, the [Nil] 205 //m 
emission could be a more convenient SFR indicator in 
the sense that it is more unambiguous, i.e. less affected 
by emissions from older stars (e.g. in the Milkyway the 
[N 11] 205 /xm flux is lower than that predicted by a linear 
L[Nii]-iiR relation when IR intensity going lower; see Fig. 
4 in Bennett ct al. 1994), and inter-comparable since the 
techniques for obtaining Ljr for low and high redshift 
objects arc so instrument-dependent. The majority of 
massive galaxies up to 2; = 5 — 6 arc already chemically 
evolved systems based on their CO and dust emissions 
(Solomon & Vanden Bout 2005) or optical/IR line spec- 
troscopy (Hamann & Ferland 1999). Therefore, Eq. (2) 
can be reasonably applicable to high- 2; starburst galaxies 
despite that the LjNiij-iiR relation possibly depends on 
metallicity. 

3.2. The scatter in the iiR-i[Nii] relation 

In Fig. 2 we can also see an indication of a steeper rela- 
tion between i[Nii] and Ljr for Ljr > 10^^ L© galaxies. 
This is more clearly demonstrated in Fig. 3 a, in which 
we plot i[Kiii]/iiR vs LiR for the same set of galaxies. 
For ULIRGs, the mean i[Nii]/-^iR (—4.69 dex) is consid- 
erably lower than that (—4.14 dex) of the galaxies with 
lower IR luminosities. Due to the small ULIRG sample 
size and the prevalence of AGN we can reach no firm 
conclusion as to whether the LjNiij/iiR decreases with 
LiR > 10^^ Lq for starburst-dominated galaxies. 

However, the scatter in the iiR-i[Nii] relation is as 
large as ~ 0.3 dex, and seems to increase with Ljr, 
specially when ULIRGs are taken into account. The 
scatter could have several different origins: metallicity 
variations and most likely, differing ionizing conditions. 
These latter effects include: the stellar UV input via the 
initial mass function, the aging of starbursts, the H 11 
rcgion/ISM geometry or the presence of an AGN. Our 
sample of galaxies contains few low metallicity objects, 
and therefore we only address ionization effects here. 

We note that the ionization parameter (U) can be es- 
pecially high in (U)LIRGs (e.g. Farrah et al. 2007, Petric 
et al. 2011), where U is defined as the dimensionless ra- 
tio of the number density of incident ionizing photons to 
the number density of hydrogen nuclei. Clearly, U can 
be measured through ionization parameter-sensitive line 
ratios, such as [O ill] /[On] and [N ill] /[Nil]. Moreover, 
theoretical models have also shown that U can be roughly 
indicated by the IR color, /eo/ZioOi by the means that 
\ogU increases with /eo/Zioo (see e.g. Abel et al. 2009). 

In Fig. 3b we plot the dependence of i[Nii]/-^iR on 
/eo/Zioo- We can see from the figure that there only ex- 
ists a very weak dependence of i[Nii]/iiR on /eo/Jioo 
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Figure 2. [Nil] 205 /jm plotted against total infrared luminosities for normal and luminous infrared galaxies. Circles (red) are galaxies 
having Herschel observations. Diamonds (blue) are galaxies from Brauher et al. (2008), for which an additional uncertainty of ~ 0.3 dex 
should be taken into account besides the plotted error bar. The superposed line represents a geometrical mean, least-squares linear fit to 
all data points except for those solid symbols. 



when log (/eo/Zioo) < -0.2. However, L[Nii]/iiR clearly 
decreases as /eo/Zioo increases when log (/eo/Zioo) ^ 
—0.2, except for a few galaxies with the warmest FIR 
colors in the sample. Moreover, the scatter in this cor- 
relation is less than that in the L[Nii]/iiR vs. Lir plot 
(see Fig. 3a) and depends little on /eo/Zioo- But we 
note that the scatter in the i[Nii]/-^iR-/6o//ioo plot is 
still large, and we hypothesize that it is caused by the 
variations in the radiation field, since N+ needs high en- 
ergy (29.6 eV) to form N"*""*", and the population of ionic 
N is determined mainly by the slope of the background 
UV spectrum, i.e. the hardness, at a given U. 

To try to account for the effects of the hardness of 
ionizing radiation field on L[nii]/Lir, we use the ratio 
of the [O III] 88 ^m emission, rather than [N III] 57 fim 
emission where there are few detections, to the [N ii] 205 
/im emission to trace the hardness of the ionizing pho- 
tons. It takes 35 eV to form which is only ~ 5 eV 
higher than that needed for the formation of N++, and 
thus the [O III] /[Nil] line ratios are even stronger indica- 
tor of the hardness. Furthermore, this ratio is insensitive 
to the gas density (Rubin 1984). Except for M82 (Duffy 
et al. 1987) and Mrk 231 (Fischer et al. 2010), all of 
the [O III] 88 fim fluxes are adopted from Brauher et al. 
(2008). 

Fig. 3c shows the relation between ./j[oiii]/-^[nii] a-nd 
-^[Nii]/-^iR- Clearly, there exists a good correlation such 
as i[Nii]/iiR decreases as i[oiii]/i[Nii] increases. There 
is only one outlier, the ULIRG Mrk 231, which shows 



both small i[iMii]/iiR and i[oiii]/i[Nii]- Excluding Mrk 
231, and using an ordinary least-squares linear fit to 
the data, we obtain log(L[Nii]/iiR) = (—3.63 ± 0.06) + 

(—0.61 ± 0.06) log(L[oiii]/i[Nii])- The Spearman rank 
correlation coefHcient p of this correlation is —0.93 at a 
> 5(7 level of significance. 

Galaxies plotted in Fig. 3 c (hereafter [Oiii] sample) 
arc also marked with red circles in Figs. 3a and 3b, 
which indicates that the hardness variation of the heat- 
ing radiation field among galaxies can largely account for 
the scatter shown in the iiR-i[Nii] relation, at least for 
LiR < 10^^ Lq galaxies. Actually, this is also supported 
by the fact that the scatter of less luminous galaxies (0.2 
dex for Brauher et al. sample; thus probably less active 
in star formation and lower U) is smaller than that of our 
luminous sample (0.28 dex). Keeping in mind that the 
conversion from [N ii] 122 /im flux to [N ii] 205 iim flux 
may contribute largely to the scatter of Brauher et al. 
sample, the real scatter in the i[Nii]"-^iR relation might 
probably be even smaller for these galaxies. 

To obtain a direct view of how much the scatter 
can be reduced, we also display the hardness-corrected 
L[Nji]/XiR, which was calculated according to the fitted 

L[Oiii]/-£'[Nii]--^[Nii]/iiR relation, for the [Oiii] sample in 
Fig. 3d (arbitrarily scaled). Comparing to Fig. 3a, one 
can see that the scatter is reduced significantly (from 
0.26 dex to 0.11 dex by excluding Mrk 231, which is still 
an outlier due to that the hardness correction has little 
effect on it, and is outside the plotted range). 
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Figure 3. Ratio of [Nil] 205 fim to IR luminosity is plotted against (a) and (d): IR luminosity, (b): IRAS 60 ^m/100 ^m color, and (c): 
[Olll] 88 /im/[Nn] 205 fim flux ratio. Symbols are the same as Fig. 2 except that here larger circles represent ULIRGs. Symbols enclosed 
by an additional red circle in panels (a) and (6) are the same sources as shown in panels (c) and (d). Eq. (1) is shown by the dashed line 
in panel (a). The ordinate of panel (d) is a nominal value after the correction for the hardness effect (using the equation shown in this 
panel). Mrk 231 is outside the plotted range of panel (d). 



All of the seven ULIRGs in our sample show some 
deficit in Ljnu] relative to Lib,. Four of the seven, in- 
cluding Arp 220, and Mrk 231, were observed by ISO 
LWS, but no [Oiii] 88 fim emission was detected at the 
3(7 limit (e.g. Brauher et a. 2008), despite all showing 
some AGN characteristics. These galaxies also show very 
warm FIR colors (see Fig. 36). Therefore, very dusty 
H II regions with high U might be needed to explain the 
deficit of [Nil] emission in these galaxies (e.g. Fischer et 
al. 1999; Abel et al. 2009; Draine 2011). 

4. SUMMARY 

In this Letter we have presented our initial results of 
the [N ii] 205 fim emission for a sample of 70 (U)LIRGs 
observed by the Herschel SPIRE FTS. Combining with 
the ISO and IRAS data, we have studied the correlation 
between L[nii] and .Z^ir, and investigated the dependence 
of L[Nii]/iiR on LiR, FIR color /eo/Aoo and the UV 
hardness indicator [O ill] /[Nil]. The main conclusions of 
our work are: 

1. For star- forming galaxies with Ljr < 10^^ Lq, 
there exists a strong correlation between L[nii] and 
LiR, namely logLjR = (4.51 ± 0.32) + (0.95 ± 
0.05) log L[Nn]: implying that the [Nil] emission 
can be used as an effective measure on SFR, partic- 
ularly for sources where iiR is difficult to measure 
such as high-z galaxies. 

2. For galaxies with ijR < 10^^ Lq in our sample, 
the L[Nii]/-^iR ratio is in the range of ~ —4.8 dex 



and ~ —3.6 dex, and varies systematically with 
the hardness/ionization parameter of the radiation 
field. However, additional effects are needed to ex- 
plain the scatter in the L^^u^-Liyi relation and their 
relatively lower L[nii]/^ir ratio for more IR lumi- 
nous galaxies in our sample. 
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